Spin properties of single electron states in coupled quantum dots 
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Spin properties of single electron states in laterally coupled quantum dots in the presence of a 
perpendicular magnetic field are studied by exact numerical diagonalization. Dresselhaus (linear 
and cubic) and Bychkov-Rashba spin-orbit couplings are included in a realistic model of confined 
dots based on GaAs. Group theoretical classification of quantum states with and without spin 
orbit coupling is provided. Spin-orbit effects on the g-factor are rather weak. It is shown that 
the frequency of coherent oscillations (tunneling amplitude) in coupled dots is largely unaffected 
by spin-orbit effects due to symmetry requirements. The leading contributions to the frequency 
involves the cubic term of the Dresselhaus coupling. Spin-orbit coupling in the presence of magnetic 
field leads to a spin-dependent tunneling amplitude, and thus to the possibility of spin to charge 
conversion, namely spatial separation of spin by coherent oscillations in a uniform magnetic field. 
It is also shown that spin hot spots exist in coupled GaAs dots already at moderate magnetic fields, 
and that spin hot spots at zero magnetic field are due to the cubic Dresselhaus term only. 

PACS numbers: 71.70.Ej,73.21.La 
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The possibility of tuning spin-orbit couplin, 
low-dimensional semiconductor electronic structures has 
stirred great interest in spin properties of lateral semicon- 
ductor electron systems in the presence of Dresselhaus 4 
and Bychkov-Rashba^ spin-orbit couplings. The for- 
mer appears in low-dimensional systems lacking inver- 
sion symmetry in the bulk (such as zinc-blende semicon- 
ductors), the latter in low-dimensional structures with 
asymmetric confining potentials. The principal question 
is what spin and charge properties and to what degree 
can be affected and manipulated by this tuning. Such 
questions are of fundamental importance for spintronics.? 

Electron spins in coupled quantum dot systems have 
been proposed to perform universal gating of quantum 
computers A The spin acts as a qubit and exchange 
coupling provides the physical realization of two-qubit 
gatesi2ii£ Another application of a controlled coupling 
between dots is spin entanglement distillation in which 
singlet and triplet states get spatially separated during 
adiabatic passage through trapped states^ Understand- 
ing of spin-orbit properties of coupled dots is thus of great 
interest to quantum information processing. 

Spin-orbit coupling provides a way for orbital de- 
grees of freedom to influence spin states. As a re- 
sult the spin dynamics is affected, making spin qubit 
operations more complex (it was shown, though, that 
two qubit operations can be performed reliably even 
in the presence of spin-orbit interaction which leads 
to anisotropic exchang o 12 ! 13 ). Furthermore, spin-orbit 
coupling leads to spin decoherence and relaxation due 
to phononsji^iSiiSiiiiS* 1 ^ limiting the operation time. 
The impressive experimental progress in coherent os- 
cillations in coupled dot systemsSSiSi* 2 ^ 3 ., as well in 
spin dephasing and spin manipulation in singl o 24 ! 25 and 
double dots2&, provides additional strong impetus for 
investigating spin states in double dots. Spin-orbit 



properties of single dots have been already extensively 
investigatedi 2 ^ 2 ^* 2 ^ 3 ^* 3 ^ 3 ^* 3 ^ 4 ^^ 3 ^ 3 ^ 3 ^ 

In this paper we investigate the role of spin-orbit cou- 
pling, represented by the Dresselhaus (both linear and cu- 
bic) and Bychkov-Rashba terms, in spin and charge prop- 
erties of two laterally coupled quantum dots based on 
GaAs materials parameters. We perform numerically ex- 
act calculations of the energy spectrum using the method 
of finite differences. We first study the general structure 
of the energy spectrum and the spin character of the 
states of the double dot system. We construct the group 
theoretical correlation diagram for the single and double 
dot states and indicate the possible transitions due to 
spin-orbit coupling. This group theoretical classification 
is used in combination with Lowdin perturbation theory 
to explain analytically our numerical results. In particu- 
lar, we show that while allowed by symmetry, the specific 
forms of the linear spin-orbit interactions do not lead to 
spin hot spots in the absence of magnetic field (spin hot 
spots are nominally degenerate states lifted by spin-orbit 
coupling 39 ). Only the cubic Dresselhaus term gives spin 
hot spots. If identified experimentally, the strength of 
the cubic term can be detected. 

We next focus on two important measurable param- 
eters: electronic g-factor and tunneling amplitude. In 
single dots the variation of the effective g-factor with the 
strength of the spin-orbit interaction has been investi- 
gated earlieriS 9 . The effect is not large, amounting to a 
fraction of a percent. Similar behavior is found for dou- 
ble dots. In our case of GaAs the contribution to the 
g-factor from spin-orbit coupling is typically about 1%, 
due to the linear Dresselhaus term. 

More exciting is the prospect of influencing coherent 
tunneling oscillations between the dots by modulating 
the spin-orbit coupling strength. Two effects can ap- 
pear: (i) the tunneling amplitude or frequency can be 
modulated by spin-orbit coupling and, (ii) the tunnel- 
ing amplitude can be spin dependent. We show how a 
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naive application of perturbation theory leads to a mis- 
leading result that (i) is present in the second order in 
linear spin-orbit coupling strengths, giving rise to an ef- 
fective tunneling Hamiltonian involving spin-flip tunnel- 
ing at zero magnetic field. Both numerical calculations 
and an analytical argument, presented here, show that 
this is incorrect and that there is no correction to the 
tunneling Hamiltonian in the second order of linear spin- 
orbit terms. The dominant correction in the second order 
comes from the interference of linear and cubic Dressel- 
haus terms. We propose to use this criterion, that the 
corrections to linear terms vanish in the second order, to 
distinguish between single and double dots as far as spin 
properties of the states are concerned. Indeed, at very 
small and very large intradot couplings the states have a 
single dot character and the correction to energy due to 
linear spin orbit terms depends on the interdot distance 
(except for the two lowest states which provide tunnel- 
ing). We find that dots are "coupled" up to the interdot 
distance of about five single-dot confinement lengths. 

In the presence of magnetic field the time reversal sym- 
metry is broken. The presence of spin-orbit coupling then 
in general leads to a spin dependent tunneling amplitude. 
Spin up and spin down states will oscillate between the 
two dots with different frequencies (for our GaAs dots 
the relative difference of the frequencies is at the order 
of 0.1%, but is higher in materials with larger spin-orbit 
coupling). This leads to a curious physical effect, namely, 
that of a spatial separation of different spin species. In- 
deed, starting with an electron localized on one dot, with 
a spin polarized in the plane (that is, a superposition of 
up and down spins), after a sequence of coherent oscil- 
lations the electron state is a superposition of spin up 
localized on one, and spin down localized on the other 
dot. A single charge measurement on one dot collapses 
the wave function to the corresponding spin state, real- 
izing a spin to charge conversion. We construct an effec- 
tive, four state (two spin and two sites) tunneling Hamil- 
tonian for the single electron double dot system, which 
takes into effect the above results. Such a Hamiltonian 
should be useful for constructing realistic model theories 
of spin dephasing, spin tunneling, and kinetic exchange 
coupling in coupled quantum dot systems. 

The paper is organized as follows. In Sec. II we in- 
troduce the model, the numerical technique, and mate- 
rials and system parameters. In Sec. Ill we review the 
benchmark case of single dots with spin-orbit coupling 
and magnetic field. Coupled double dots are studied in 
Sec. IV, separately in zero and finite magnetic fields. We 
conclude with the discussion of our results in Sec. V. 



II. MODEL 

We consider conduction electrons confined in a [001] 
plane of a zinc-blende semiconductor heterostructure, 
with additional confinement into lateral dots given by 
appropriately shaped top gates. A magnetic field B is 



applied perpendicular to the plane. In the effective mass 
approximation the single-electron Hamiltonian of such a 
system, taking into account spin-orbit coupling, can be 
decomposed into several terms: 

H = T + Vc + H z + H B R + H D + H D3 . (1) 

Here T = ti 2 K 2 /2m is the kinetic energy with the ef- 
fective electron mass m and kinetic momentum KK. = 
Kk + eA = —ihV + eA; e is the proton charge and A = 
B(-y /2,x/2,0) is the vector potential of B = (0,0,5). 
Operators of angular momentum with mechanical and 
canonical momenta are denoted as L = r x (KK) and 
1 = r x (tik). The quantum dot geometry is described 
by the confining potential Vc(r). Single dots will be de- 
scribed here by a parabolic potential Vc = (l/2)mu) 2 l r 2 , 
characterized by confinement energy Eq = Kljq/2 and 
confinement length 1$ = (ft/mwo) 1 ^ 2 , setting the energy 
and length scales, respectively. Coupled double dots will 
be described by two displaced (along x) parabolas: 

V$ d = ±mu 2 l(\x\-l df+y 2 }; (2) 

the distance between the minima is 2d measured in the 
units of Iq. The Zeeman energy is given by Hz = 
(g* /2)fisO'zB, where g* is the conduction band g- factor, 
fis is the Bohr magneton, and a z is the Pauli matrix. In 
order to relate the magnetic moment of electrons to their 
orbital momentum, we will use dimensionless parameter 
az = g*m/2m e , where m e is the free electron mass. 

Spin-orbit coupling gives additional terms in confined 
systems^ The Bychkov-Rashba Hamiltonian, 5-6 

Hbr = otBR (cr x Ky — a y K x ) , (3) 

appears if the confinement is not symmetric in the growth 
direction (here z). The strength olbr of the interac- 
tion can be tuned by modulating the asymmetry by top 
gates. Due to the lack of spatial inversion symmetry in 
zinc-blende semiconductors, the spin-orbit interaction of 
conduction electrons takes the form of the Dresselhaus 
Hamiltonian 4 which, when quantized in the growth direc- 
tion z of our heterostructure gives two terms, one linear 
and one cubic in momentum^ 

H D = lc (K 2 )(-a x K x + a y K y ), (4) 
H D3 = ( 7c /2) (a x K x K 2 - a y K v Kl) + He, (5) 

where 7 C is a material parameter. The angular brackets 
in Hr, denote quantum averaging in the z direction-the 
magnitude of Hp depends on the confinement strength. 
We will denote the sum of the two linear spin-orbit terms 
by Hn n = Hb> + Hbr- The complete spin-orbit cou- 
pling is then Hso — Hin + Hr>3- We find it useful 
to introduce dimensionless strengths of the individual 
terms of the spin-orbit interaction by relating them to 
the confinement energy of a single dot Eq. We denote 
&br. = obr/EqIq anda£> = ^ydk 2 } / EqIq for linear terms, 
and ao3 — 7c/2£'o'o f° r the cubic Dresselhaus term. 
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In our numerical examples we use Eq = Huiq/2 = 1.43 
meV for the confinement energy, which corresponds to 
the confinement length of Iq = 20 nm. We further 
use bulk GaAs materials parameters: m = 0.067 m e , 
g* = -0.44, and 7 C = 27.5 eVA 3 . For (k 2 ) we choose 
5. 3x 10 _4 A 2 , which corresponds to 7c(fcf) = 14.6 meV A. 
This value of (k 2 ) corresponds to the ground state of a 6 
nm thick triangular potential well^S For a br we choose a 
generic value of 4.4 meVA, which is in line of experimen- 
tal observations The dimensionless parameter of the 
Zeeman splitting then is az = —0.015, while the relative 
strengths of the spin-orbit interactions are olbr ~ 0.015, 
a_D ~ 0.05, and ao3 ~ 0.001. Except for anti-crossings, 
the spin-orbit interaction is a small perturbation to the 
electronic structure; it is, however, essential for investi- 
gating spin properties. 

Our analytical calculations will often refer to the Fock- 
Darwin^a* spectrum, which is the spectrum of Hamilto- 
nian for a single dot with Hso = 0. The correspond- 
ing wave functions \& (expressed in polar coordinates r 
and 0), and energies e are 

* n ,iAr,4>)=CpMe~» 2 / 2 LW(p 2 y l H(«), (6) 

I 2 he 
e nM =2E -l-(2n + \l\ + 1) + B—(l + a z a),(7) 
l B 2m 

where p — t/Ib is the radius in the units of the effec- 
tive confinement length Ib, defined by l 2 B — ^/V(l + 
B 2 e 2 l\i '4h 2 ); n and I are the principal and orbital quan- 
tum numbers, respectively, C is the state dependent nor- 
malization constant, and l)n are associated Laguerre 
polynomials. Spinors £(cr) describe the spin a state of 
the electrons. Since the parabolic dot has rotational sym- 
metry in the plane, the states have well defined orbital 
momentum I and spin a in the z-direction. We also in- 
troduce a useful dimensionless measure 9 of the strength 
of the magnetic field induced confinement compared to 
the potential confinement: 9 = Bel 2 B /2h, < 9 < 1. The 
parameter 9 gives the number of magnetic flux quanta 
through a circle with radius I b ■ For large magnetic fields 
(Bel 2 /2h) 2 /2. The confining length can be ex- 
pressed as Ib = ^o(l — 9 2 ) 1 ^ 4 . 

As it is not possible to solve for the spectrum of 
Hamiltonian analytically, we treat it numerically with 
the finite differences method using Dirichlet boundary 
conditions (vanishing of the wave function at bound- 
aries). The magnetic field is included via the Peierls 
phase: if H(Yi,Yj) is the discretized Hamiltonian con- 
necting grid points and Yj at B = 0, the effects 
of the field are obtained by adding a gauge phase: 
H(ri,Yj) exp[i(e/h) J 3 A.dl]. In our simulations we typ- 
ically use 50 x 50 grid points. The resulting matrix eigen- 
value problem is solved by Lanczos diagonalization. The 
achieved accuracy is about 10~ 5 . 



III. SINGLE DOTS 

As a starting point we review the effects of spin-orbit 
coupling in single dots. We are interested in the changes 
to the spectrum and, in particular, to the magnetic mo- 
ment of the lowest states, that is, to the effective g-factor. 
The calculated spectrum of a single dot is shown in Fig. 
^ There are three ways in which spin-orbit coupling 
influences the spectrum: (i) First, the levels are uni- 
formly shifted down, in proportion to a 2 (by a here we 
mean any of cxbr, ctp, or olds)- (h) Second, the de- 
generacy at B — is lifted, again in proportion to a 2 
^p). (iii) Finally, at some magnetic field the level cross- 
ing of the Fock-Darwin levels is lifted by spin-orbit cou- 
pling. The resulting level repulsion is linear in a Q;). 
The states here are the spin hot spots, that is states in 
which both Pauli spin up and down species contribute 
significantlyi 19 i 32 i 39 

The above picture can be understood from general 
symmetry considerations within the framework of pertur- 
bation theory. All spin-orbit terms commute, at B = 0, 
with the time inversion operator T — ia y C , where C is 
the operator of complex conjugation. Therefore Kramer's 
degeneracy is preserved so that the states are always dou- 
bly degenerate. The linear terms can be transformed into 
each other by a unitary transformation (cr x + a y ")j\[2 
(spin rotation around [110] by 7r ), which commutes with 
Hq ■ Therefore the effects on the energy spectrum induced 
individually by the linear Dresselhaus and the Bychkov- 
Rashba terms are identical at B = Q. At finite magnetic 
fields the two interactions give qualitatively different ef- 
fects in the spectrum, especially for spin hot spots, as 
discussed below. 

For any B the following commutation relations hold 
for the linear terms: 

[H B R,lz + s z ] =0, [H D ,l z -s z ]=0. (8) 

This means that Hbr commutes with the angular mo- 
mentum, while Hjj does not. This will influence the in- 
terference between the two terms in the energy spectrum. 
We can use the Fock-Darwin states as a basis for pertur- 
bation theory. Up to the second order the energy of state 

\i) = ^n,l,a IS 

jp _l u\u u\ i V" (i\ H so\j)(j\Hso\i) , Q . 
Ej = €j + {i\H so \i} + } j . (9) 

The first order correction is zero for all spin-orbit terms 
since H$o contain only odd powers of K whose expec- 
tation values in the Fock-Darwin states vanish. If the 
perturbation expansion is appropriate, the spin-orbit in- 
teractions have a second order (in a) effect on energy. 

Both linear spin-orbit terms couple states with orbital 
momenta I differing by one. It then follows from the 
commutation relations JSJ that Hbr preserves the to- 
tal angular momentum I + s, while Hp. preserves the 
quantity I — s. As a result, there is no correction to 
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FIG. 1: Energy spectrum of a single dot in magnetic field, a) 
The Fock-Darwin spectrum, Eq. J7[. b) Lowest orbital ex- 
cited levels (n = 0, |i| = I) without (dashed) and with (solid) 
spin-orbit coupling. Arrows indicate the spin states. For clar- 
ity the energy's origin here is shifted relative to case a). Both 
the shift in energy levels as well as the splitting at B = 
grow as a 2 , c) Anti-crossing at the critical magnetic field 
(here about 13 T). For clarity, a linear trend was subtracted 
from the data. 



the energy from the interference terms between Hbr 
and H D in Eq. ©: {i\H BR \j)(j\H D \i) = 0. As for 
the cubic Dresselhaus term, only the following orbital 
states are coupled: (I, |) — > {(I + 3,i),(Z — 1,1)} and 
(hi) — * {(!■ ~ 3, f), (Z + l,t)}- Due to these selection 
rules there are no interference terms ~ Hd^Hbr, but 
terms ~ Hb?,Hb will contribute to energy perturbation. 
The Bychkov-Rashba and Dresselhaus Hamiltonians act 
independently on the Fock-Darwin spectrum (up to the 
second order). 

To gain more insight into the perturbed structure of 
the spectrum at B = 0, we rewrite Eq. 10 using an aux- 
iliary anti-hermitian operator H^ defined by the com- 
mutation relation [Hq, Hg P ] — Hso- If such an operator 



E 



(i\H so \j)(j\Hso\i) 



= (i\~[HZ,Hso}\i) + 



-Re{(i\H S oPNH° s p \i)), 



(10) 



where Pv is the projector on the subspace M of the states 
excluded from the summation. In our case here it is just 
one state, AT = The last term in (|10fl then vanishes. 

The auxiliary operator for H03 is not known and if found, 
it must depend on the confining potential. Operators for 
the linear terms are 4^ 



D 

n BR 



-i(a D /2l )(xa x - ya y ), (11) 
i(a BR /2l )(ya x - xoy). (12) 



The corresponding commutators are (in the zero mag- 
netic field K = k, L z — l z , 8 — 0; the last expression will 
be useful later) 



[H°v,H D ] = -E al(l-a z L z ), 
[H° B P R ,H BR ] - 



(13) 
(14) 



-E a 2 BR (l + a z L z ), 
H D P ,H D3 ] = E ll aD a m (Kl + K 2 y - 

-2o z \xK y Kl - yK x Ky - 2i6(xK x + yK y )]) .(15) 



Because [Hp \ Hbr] + [H° b p Ri Hd\ = 0, the corrections 
to the second order perturbation add independently for 
Hbr and Ho (as also noted above from the selection 
rules), we can introduce H°^ — H° D V + H B P R . An alter- 
native route to Eq. (|10f> is to transform the Hamiltonian 
within U = exp(-H° p ) toH = H - (l/2)\H S0 ,Hf ] 
in the second order of a. The final result can be also 
obtained in a straightforward way by using the Thomas- 
Reiche-Kuhn sum rule in the second order of perturba- 
tion theory with the original spin-orbit terms. The re- 
sulting effective Hamiltonian is (terms depending on ckd3 
are omitted here) 

H = Ho-E {a 2 D +a 2 BR )/2+E a z L z (a 2 D -a% R )/2. (16) 

This Hamiltonian, in which the spin-orbit coupling ap- 
pears in its standard form, neatly explains point (ii) 
about the lifting of the degeneracy at B = 0. The levels 
in Fig. ^h for example, are four fold degenerate (|/| = 1, 
|cr| = 1) without spin-orbit coupling. Turning on, say, 
Hb, will split the levels into two groups: energy of the 
states with la > would not change in the second order, 
while the states with la < will go down in energy by 
E n a D , as seen in Fig. ^p. 



A. Spin hot spots 

Spin hot spots are states formed by two or more states 
whose energies in the absence of spin-orbit coupling are 
degenerate or close to being degenerate, while turning on 
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FIG. 2: Calculated corrections to the effective g-factor by 
spin-orbit interactions. Formulas 11811 scaled by the values 
at B = (and thus independent on aso) are plotted. The 
actual numerical values of Sg at B — are: SgD-D(0) = 
l.OxHT 2 , Sg B R-B R (0) = 8.6xHT 4 , Sg D - D3 (0) = 9.4xl0" 4 , 
5go3-D3(0) — 2.5xl0~ 5 . At the anti-crossing 5gD_D(Bacr) = 
2.4 x KT 3 , SggR-BRiB^) = 1.0 x HT 4 , Sgo-D^B^) = 

3.4 x 10" 4 . 
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the coupling removes the degeneracy^ Such states are 
of great importance for spin relaxation, which is strongly 
enhanced by their presence The reason is that the 
degeneracy lifting mixes spin up and spin down states 
and so transitions between states of opposite magnetic 
moment will involve spin flips with a much enhanced 
probability compared to normal states. 

Figure shows an interesting situation where two de- 
generate levels are lifted by spin-orbit coupling>i2i22 The 
lifting is of the first order in a, unlike the lifting of de- 
generacy at B = in which case the degenerate states 
are not directly coupled by H$o- In a finite magnetic 
field, at a certain value B a cr, the states of opposite spins 
and orbital momenta differing by one cross each other, 
as follows from the equation J7J). The crossing field is 
Sacr ~ \a z \' 1/2 h/{ell), which is about 13.4 T for our 
parameters (making the confinement length larger the 
magnitude of the field would decrease). Spin-orbit inter- 
action couples the two states thereby lifting the degener- 
acy. For GaAs, where g* < 0, only the Bychkov-Rashba 
term couples the two states. The Dresselhaus terms are 
not effective {Hb>3 would introduce such a splitting at 
3i?acr)- The energy splitting due to Hbr is 



A = cE OL B r{1 — #acr)(l — #acr) 



,1/4 



(17) 



where c, which is a number of order 1, depends on the 
quantum numbers of the two states. Since spin hot spots 
at -Bacr are due only Hbr, the splittings could help to 
sort out the Bychkov-Rashba versus Dresselhaus contri- 



butions. Figure HI 1 shows the calculated level repulsion 
for states n = 0, 1 = 0, a ={ and n = 0,1 = —1,(7 =|. 
The magnitude of A, though being linear in czbr, is on 
the order of 10 -3 meV and thus comparable to the energy 
scales associated with quadratic spin-orbit perturbations. 



B. Effective g-factor 

When probing spin states in quantum dots with mag- 
netic field, important information comes from the mea- 
sured Zeeman splitting. We will focus here on the two 
lowest spin states and calculate the effective g-factor as 
g = (-B ,o,| - E qm )/{hbB). If H S o = 0, then the effec- 
tive g-factor equals to the conduction band value g*. The 
actual value in the presence of spin-orbit coupling is im- 
portant for understanding single spin precession in mag- 
netic field, which seems necessary to perform single qubit 
operations in quantum dots. We have obtained the fol- 
lowing contributions to the g-factor from non-degenerate 
(that is, excluding spin hot spots) second-order pertur- 
bation theory [Eq. ©] (for linear spin-orbit terms these 
are derived also in 29 ' 4 -): 
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8m/m e 6Vl - 9 2 \ l-9(l+2a z ) 
(1 - 6>) 4 (1 + Of -3(1 - ef 
3-9(l + 2a z ) + 3-9(3-2a z ) + 

3(1 + 6Q 6 2(l + 6>) 2 (l + 6> 2 ) 2 

3 + 6>(3 - 2a z ) 1 + 0(1 + 2a z ) 

(1-0) 2 (1 + 0) 4N 



[l + 2a z ) J' 



(18) 



Here SgA-B stands for a correction that is proportional 
to olaolb- 

The functions (|18fl are plotted in Fig. [21 We can 

understand the limits of 5g at B — > oo (9 — > 1) if we 
notice that in the natural length unit Is the momen- 
tum K x = -id x - yBe/2h = ^ 1 [-i9 x /i B - 0{y/l B )]- In 
the limit B — > oo the matrix elements of Hd, which is 
linear in K, scale as lg 1 , while the Fock-Darwin ener- 
gies scale as lg 2 . The second order D-D correction to 
©0,0,4, — -©o.o, T is thus independent of Is', it converges to 
— E a 2 D /(l + a z ). The BR-BR correction is analogous, 
with the limit E a 2 BR /(l — a z ). To get the g-factor we di- 
vide the energy differences by hbB and get figD-D (9 —> 
1) oc B~ x ; similarly for H B r- Since H D3 scales as lg 3 
one gets SgD-D3 {9 — > 1) — > 2a£>OD3m/(l + a z )m e and 
figD3-D3(9 — > 1) oc B, It seems that for increasing B 
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magnetic field 


SO terms 


symmetries of H 
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under I x , I y 


numbers for 


n.l.cr 

9i 


5 = 


none 






sentation 


transforms 


I - even 


/ - 


odd 




BR 


—ia x I x ,—ia y Iy,—ia z I,T 






as 


L D 


L D 




D,D3 


—ia y I x ,—icr x Iy,—icr z I,T 




Ti 


1 


1 1 




-1 -1 




all 


—ia z I,T 




r 2 


X 


-1 -1 




1 1 


B > 


none 


—io z I,R z 




r 3 




-1 1 




1 -1 




any 


—io z I 




r 4 


y 


1 -1 




-1 1 



TABLE I: Symmetries of the double dot Hamiltonian for dif- 
ferent spin-orbit terms present at B = and B > 0. Here 



Ix(Iy) means reflection x 



-x (y 



-y), I 



R z = exp(— i(f>a z /2) is the rotation of a spinor by angle <j> 
around z-axis; R n is a spinor rotation around an arbitrary 
axis n and T is the time reversal symmetry. The identity 
operation is not listed. 



there inevitably comes a point where the influence of 
Hd3 on the g-factor dominates. But at B = £? a cr there 
is an anti-crossing of the states (0,0, J.) and (0, — 1, f) 
so for larger B the g-factor does not describe the en- 
ergy difference between the two lowest states, but be- 
tween the second excited state and the ground state. 
The value of B where <5gr>3-.D3 = ^9d-d is given by 
B w {hj el\)au I OL^\pl. For GaAs parameters it is ~ 25 
T. After anti-crossing the first exited state is ^o.-id arL d 
its energy difference to the ground state is divided by 
HbB is oc 1/B 2 for Hq. The corrections from spin-orbit 
terms are D-D, BR-BR oc 1/B 5 , D-D3 oc 1/B 4 , and 
D3-D3oc 1/B 3 . 



IV. DOUBLE DOTS 

A double dot structure comprises two single dots close 
enough for their mutual interaction to play an important 
role. Here we consider symmetric dots modeled by 
of Eq. (J2J. Such a potential has an advantage that in 
the limits of small d — * and large d — ► oo, the solutions 
converge to the single dot solutions centered at x = and 
Hod, respectively. We denote the displaced Fock-Darwin 
states as *^ jCT 0, y) = ^ n ,i,a(x ± l d, y). 

The symmetries of the double dot Hamiltonian with 
spin-orbit couplings are listed in Tab. [fl The time sym- 
metry is always present at B — 0, giving Kramer's dou- 
ble degeneracy. The rotational space symmetry from the 
single dot case is lost; instead there are two discrete sym- 
metries - reflections I x about y and I y about x. In zero 
magnetic field and without spin-orbit terms, the Hamil- 
tonian has both I x and I y symmetries. If only Rashba 
or Dresselhaus terms are present, we can still preserve 
symmetries I x and I y by properly defining the symmetry 
operators to act also on the spinors (forming the double 
group). The Bychkov-Rashba term, Hq + Hbr, is invari- 
ant to operations defined by the spatial invariance. This 
is not the case of Hp, since here the operators —ia y I x 
and —ia x I y do not describe a spatial reflection of both 
the orbital and spinor parts. The symmetry operations 



TABLE II: Notation and transformation properties of Civ 
representations. L and D are the coefficients of the depen- 
dence of g"' l,a on the single dot functions (see text). 



for Hbr and Hp are connected by the unitary transfor- 
mation (a x + (T y )/y/2, which connects the two Hamilto- 
nians themselves. Finally, if both spin-orbit terms are 
present, or at B > 0, the only space symmetry left is 
I I x I y . 

In the following we consider separately the cases of zero 
and finite magnetic fields. 



A. Energy spectrum in zero magnetic field, 
without spin-orbit terms. 

If no spin-orbit terms are present the group of our 
double dot Hamiltonian is C 2v ® SU(2). The SU(2) 
part accounts for the (double) spin degeneracy. The or- 
bital parts of the eigenstates of the Hamiltonian therefore 
transform according to the irreducible representations of 
Civ- The representations^ Ti, % = 1...4, along with their 
transformation properties under the symmetries of Ci v , 
are listed in Tab. ^ The symmetry properties will be 
used in discussing the perturbed spectrum. 

We denote the exact eigenfunctions of the double dot 
Hamiltonian with Hso — as Tf£ where a(b) is the 
single dot level to which this eigenfunction converges as 
d — > (oo); i labels the irreducible representation, a de- 
notes spin. We have chosen the confining potential to 
be such, that at d — ► 0(oo) the solutions of the double 
dot Ho converge to the (shifted) Fock-Darwin functions, 
if properly symmetrized according to the representations 
of Civ These symmetrized functions will be denoted as 



where (up to normalization) 



The numbers Dj(Lj) for different irreducible representa- 
tions are in the Tab. [H] 

Generally, up to normalization, the exact solution can 
be written as a linear combination of any complete set of 
functions (we omit the spin index which is the same for 
all terms in the equation) 



n ,l 



+ c ( n ' 03* 



(20) 



The last equation indicates the fact, that for a function 
r° b in the limit d — > 0(oo), there will be a dominant g- 
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function in the sum with the numbers no, Iq given by the 
level a(b) and the coefficients c for the other functions will 
converge to zero. We term the approximation c(n, I) — 
as a linear combination of single dot orbitals (LCSDO). 

Knowing the representations of the double dot Hamil- 
tonian and the fact that Fock-Darwin functions form 
50(2) representations (reflecting the symmetry of single 
dot Hq) we can decompose all single dot levels into the 
double dot representations and thus formally construct 
the energy spectrum of a double dot using the symmetry 
considerations only. 

Following the standard technique for constructing such 
correlation diagrams (connecting states of the same rep- 
resentation and avoiding crossing of lines of the same rep- 
resentation) we arrive at the spectrum shown in Fig. [3] 
The ground state transforms by the symmetry operations 
according to ri (identity), while the first excited state ac- 
cording to T2 (x). This is expected for the symmetric and 
antisymmetric states formed by single dot ground states. 
The symmetry structure of the higher excited states is 
important to understand spin-orbit coupling effects. In- 
deed, the spin-orbit terms couple two opposite spins ac- 
cording to certain selection rules. Since Hd, for example, 
transforms similarly to iffij/, it couples the ground state 
Ti with T2 and T4. In general, odd numbered representa- 
tions can couple to even numbered representations. The 
same holds for Hbr and Hps- If we include either Hbr 
or Hd into the Hamiltonian, and consider spinors as the 
basis for a representation, the states would transform ac- 
cording to rg, the only irreducible representation of the 
double group of Civ 

The calculated numerical spectrum for our model 
structure is shown in Fig. 0] There is a nice qual- 
itative correspondence with Fig. |3| In Fig. 01 by 
vertical bars we denote coupling through Hjj or Hbr 
(\{i\H D \j)\ = \{i\H BR \j)\). The couplings follow the se- 
lection rule described above. Since there are several level 
crossings in the lowest part of the spectrum, a question 
arises if spin hot spots are formed in the presence of spin- 
orbit coupling. It turns out, that there is no first-order 
level repulsion at the crossings due to the linear spin- 
orbit terms because the levels are not coupled by the 
linear terms, even though such couplings are allowed by 
symmetry. There are no spin hot spots due to the lin- 
ear spin-orbit terms at zero magnetic field. For example 
r^ 1 and rf 1 are not coupled by spin-orbit terms, and 
therefore their degeneracy (at 2dlo « 50 nm) is not lifted 
by linear spin-orbit terms as we would expect from sym- 
metry (actually, there is an anti-crossing which is of the 
order af in , instead of the expected ou n ). The cubic Dres- 
selhaus term gives here (and also in other crossings that 
conform with the selection rules) a linear anti-crossing, 
as one expects. The absence of anti-crossings from the 
linear spin-orbit terms will be explained in the next sec- 
tion. 

Since our main goal here is to study the effects of spin- 
orbit coupling on the tunneling between the two dots, we 
first look at the tunneling for if so = 0. The tunneling 
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FIG. 3: Single electron spectrum of a symmetric (C2v) lateral 
double dot structure as a function of the interdot separation, 
at B = 0, derived by applying group theoretical considera- 
tions. Single dot states at d = and d = 00 are labeled by 
the principal (n) and orbital (I) quantum numbers, while the 
double dot states are labeled according to the four irreducible 
representations V% of Can. The lowest double dot states have 
explicitly written indices showing the excitation level of the 
d — and d = 00 states they connect. Every state is doubly 
(spin) degenerate, and spin index is not given. 



energy, SEt /2, is given by the difference between the ener- 
gies of the symmetric ground state and the asymmetric 
first excited orbital state: 5E t = Ea — E$- We com- 
pare the LCSDO approximation with our numerically 
exact calculations. The functions Eq. I|19|l are not or- 
thogonal. If we introduce the overlap integrals between 
these functions (all the indices of a g-function are denoted 
by one letter here) by Sij = (i\j), and the Hamiltonian 
matrix elements Hij = (i\H\j), variational theory gives 
for the expansion coefficients a of a double dot state 
r = J^Cjli) the generalized eigenvalue equation, 



(21) 



We will compute the energy of the two lowest orbital 



double dot states, r"" 



per 

1 s> 



plO 

1 2(7 



A' 



which are a 



symmetric and an antisymmetric state with respect to x. 
The energies are denoted as Eg , Ej± , where index zero 
indicates the absence of spin-orbit coupling. 

We first use the LCSDO approximation suitable for the 
limit d — > 00, that is we approximate F s 



0,0.0- 

92 ■ 



0.0, cr , 

9i , and 



This means the basis in Eq. (|21|l consists 



of one function and the solution for the energy is Ei 
Hii/Su. For the considered states we obtain: 



E 



(0) 



E 



(0) 



2E 



2En 



1 + [1 - 2d/v^F] 



d 2 Erfc(d) 



1 + e~ d2 
d 2 Evic(d) 



-a 2 



(22) 
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FIG. 4: Calculated energy spectrum of a double quantum dot 
at B — 0, as a function of interdot distance. Spin-dependent 
terms are not included. Vertical bars indicate couplings due 
to spin-orbit coupling. Group theoretical symbols are shown 
with the lines on the left. Single dot levels are denoted by 
the highest orbital momentum (0, 1, 2, ...) present in the 
degenerate set. This labeling is on the right. Every state is 
doubly degenerate, and spin index is not given. 
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FIG. 5: Calculated tunneling energy 8E[ 0) for the double dot 
Hamiltonian Ho- Numerical results are compared with the 
analytical (full) expression for SE^ — — where the 
energies of the asymmetric and symmetric states are given by 
Eq. 1221 . The leading order expression is that of Eq. 1231 . 
valid for large d. The inset shows the relative error (with 
respect to the numerical calculation) of the two analytical 
results. 



In the limit of large interdot separation (that is we make 
expansion in powers of e~ d and take the first term in 
this expansion as the leading order), Eg)^ w E a [2 ± 
(2/y / 7r)dexp(— d 2 )], where the minus (plus) sign is for 
E^ (E { ° ] ). The tunneling energy then is 

5E[ 0) ^E ^de- d2 . (23) 
v 7r 

To understand this result, we get it once again by tak- 
ing a two dimensional basis in Eq. (|21|l consisting of 
functions 'Fq a and 'Fq^q ct . In the limit d — > oo both di- 
agonal elements of the matrix H converge to the energy 
of the Fock-Darwin ground state (2Eo) . Then the differ- 
ence of the eigenvalues of H is given by the off-diagonal 
matrix element, which is H12 = ('Fq al^^o a) ~ 

Any further refinement beyond LCSDO would not con- 
tribute significantly to the calculated SE^°\ For ex- 
ample, for the ground state to go beyond LCSDO, we 
take the next excited function of the same symmetry 
and get the basis for Eq. (|2"T|l to be gl' 1 ^}- 
The off-diagonal element in the matrix H is H12 — 
{g^ \H \g^) fa E P(d)e- d \ where P(d) is a polynomial 
in d. Since now ff u — H 2 2 is of the order of Eq, the 
correction from non-diagonal terms will be of the order 
of (e~ d ) 2 and will not change the leading order result. 

The numerical result, together with the analytical 
forms of 5Et [one using the complete expressions Eq. 
(|22p. the other for the limiting case of large d, Eq. I|23|) ]. 
is plotted in Fig [3] The complete expression is in excel- 
lent agreement with the numerics, over the whole range 



of d, including the limit d — > 0. As for the leading or- 
der expression, it becomes an excellent description for 
the tunneling energy at distances roughly twice the dot 
confinement length (40 nm in our case) , as seen from the 
inset to Fig. [SJ 

Using a LCSDO approximation <?i'°' CT for the ground 
state is correct for both limits d — ► and d — > 00, be- 
cause r°° converges to the single dot level in both lim- 
its. Therefore it is reasonable to expect that the ap- 
proximation will be equally good for the whole range 
of d. However, for the first excited state T^ the two 
limits go into different single dot levels and the proper 
LCSDO approximations for this state are g^ 0,a and g^ l,a 
in the limit d — > 00 and d —> respectively. However 
^0,0, <j ^ g^ l,IJ 4- g ) i ,l,IJ as d — * and thus for both ground 
and first excited states using a LCSDO approximation 
suitable for d — * 00 gives good results in the whole range 
of d. We will see, that this will not be true for a non- 
zero magnetic field and describing by g 2 ' ' a will give 
a much higher error. A remedy is to go beyond LCSDO 
for r^ , for example by taking the base for Eq. I|21l) to 
include two g-functions, each correct in one of the limits 
d — > 0, d — * 00. We will not present the results of such 
computations since formulas become more complicated 
without giving better understanding. 

Higher orbital states can be treated similarly. Starting 
from level 2 there are more functions of the same repre- 
sentation in one single dot level, therefore the basis for 
the Hamiltonian Eq. (|21|l giving the leading order must 
contain more that one function. 
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FIG. 6: Calculated corrections to the tunneling energy SE t 
from spin-orbit terms at B = 0. The labels indicate which 
spin-orbit terms are involved. Only D-D3 and D3-DS are 
of second order. The remaining contributions are of fourth 
order. 



B. Corrections to energy from spin-orbit coupling 
in zero magnetic field. 

When we add Hso to H , the structure of the correc- 
tions to the energies of the two lowest states up to the 
second order in spin-orbit couplings can be expressed as 



E 



(2) 



-Ai(a 



B t a 2 



D3 



+ Cia D a D3 , (24) 



where i is either 5* or A. The coefficients A,B, and C 
are positive for all values of the interdot distance. The 
differences Aa~ As, ■ ■ ■ approach zero as d — > oo. We will 
argue below that As = A a = 1/2 with the exception of a 
very small interdot distance (less than 1 nm). There are 
thus no contributions from the linear spin-orbit couplings 
to 5E t in the second order. Only the cubic Dresselhaus 
term contributes, either by itself or in combination with 
the linear Dresselhaus term. Spin-dependent tunneling 
is greatly inhibited. 

Numerical calculation of the corrections to 5E t from 
spin-orbit couplings are shown in Fig EJ The dominant 
correction is the mixed D-D3 term, followed by D3-D3. 
These are the only second order corrections. For GaAs, 
and our model geometry, these corrections are about 4 
and 5 orders of magnitude lower than SE^ . The correc- 
tions, when only linear spin-orbit terms are present, are 
much smaller since they are of the fourth order. The dra- 
matic enhancement of the corrections from linear spin- 
orbit terms close to d = is due to the transition from 
coupled to single dots. We will explore this region in 
more detail later. 

We will first show how naive approaches to calculating 
spin-orbit contributions to tunneling fail to explain the 
above results. We use the example of the linear Dres- 
selhaus term. The simplest way to include this term 



is to begin with the two lowest orbital states (that are 
four states including spin), g{ ' CT and g 2 : • Because of 
the time reversal symmetry the resulting 4x4 matrices 
H and S from Ea. (|21|) block diagonalize into two equal 

2x2 matrices with elements Hu 
and H 12 = (5?'°' T 



E 



(0) 



Ho 



E 



(0) 



22 - ^ A 
,2 

iEoaode" ; S is the 
since the two states are orthogonal 

(0) 



\H D \gl M 



unit matrix now, 

due to symmetry. Using the large d limit for 8E. 
Eq. (|23[) . we obtain the perturbed energies E$(a) — 
2Eo±EoyJ '4/7T + a 2 D de~ d with the minus (plus) sign tor- 
s' (A) . In the second order of old the symmetric and an- 
tisymmetric level energies have opposite contributions, 
giving SE t « [(4£ /V^) + (E V^/2)a 2 D ]dexp(-d 2 ), in 
contrast to the numerical results where there is no de- 
pendence on a 2 D in the second order. Enlarging the basis 
by the first excited orbital states (all together 12 states 
including spin), that is, include ,a , the symmetric and 
antisymmetric level energies will have the correct limit for 
the spin-orbit contributions, —Eqol d J2, at d — > 00. At 
finite values of d the difference from this limit value is 
less than 2%. That means there are still terms of order 
a 2 D in SE t . Could using a renormalized basis help? We 
could, for example, use symmetrized states of the sepa- 
rated dots that include spin-orbit terms. It is not difficult 
to see that this would not work either: the perturbed 
single dot ground state, for example, contains the spin 
admixture from the first excited orbital states. This is 
then similar to using the 12 state basis in the variational 
approach. 

From the previous example one can see that to get a 
correct spin-orbit contribution to the energy of a state, 
it is not enough to include just a few terms in the sum in 
Eq. @. Instead we employ the operators H op given in 
Eqs. i|llll2|) . To get a contribution for a particular state, 
say \i), we apply the Lowdin perturbation theory^ For 
this one has to identify states \j) which are degenerate 
with |i) with respect to the perturbation Hso and which 
have to be treated exactly. The rest of the states can be 
treated perturbatively. The condition for a degeneracy 



of two states can be taken as \E, 



(0) 



E 



(o) 



< 



OtSO,Oi\h 



when one considers linear and cubic terms respectively. 
The finite set of the degenerate states will be denoted by 
N ' . The effective Hamiltonian H cS acting in Af is 



1 x 



(H a 

^(Hso 



Hso)ij + 

' k(Hso)kj 



{Hso)ik(Hso)kj 



2 ^ L 

k£N 



E 



(o) 



E 



(0) 



E 



(o) 



E. 



(o) 



(25) 



For the example of the linear Dresselhaus term, we can 
now use Eq. Ill Oil and Ijl3|l to obtain 



H c t f = (H + H D ) l3 - -a 2 D E Q (1 - + R tJ , (26) 



where 



Rij = 2 m D P^H° D p - H° D P P M H D \j) . (27) 
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First we note that existence of the operator H°p means 
that the coupling through Hp between any two states 
is always much smaller then the difference of the un- 
perturbed energies of these two states, since (Hp)ij — 
{Ef ] - E^)(H° D p ) l3 ~ {Ef ] - Ef ] )a D . Then one can 
partially diagonalize the effective Hamiltonian to elimi- 
nate the off-diagonal Hp terms. It turns out, that this 
leads to a cancellation of the terms Ho and R. The ef- 
fective Hamiltonian is then 

Htf = (flb)« - \a 2 D E (1 - v s l,).. . (28) 

This completes the way to get Eq. (|lfil) using Lowdin 
perturbation theory. There are no linear effects on the 
double dot energy spectrum from linear spin-orbit terms, 
which explains the absence of spin hot spots even though 
symmetry allows that. 

The spin-orbit interaction can influence the energy 
only through the operator l z , which is of the represen- 
tation T3, from where we get selection rule-the allowed 
coupling is between functions of representations Ti - T3 
and ^-IV Looking at Fig. 21 accidental degeneracies 
of states with such representations are not present in the 
lower part of the spectrum. The crossing of T^ 1 with T^ 1 
considered in the discussion to Fig. 0]also does not fol- 
low the selection rule, hence why the anti-crossing is of 
the third order. From the selection rules one can imme- 
diately see that also the expectation value of l z is zero 
in any state. This result is more general and holds also 
if the symmetry of the potential is lower (or none), since 
it follows from the fact that the Hamiltonian Hq is real, 
so one can choose eigenfunctions to be real. Then the 
expectation value of any imaginary operator, such as l z , 
must vanish. We conclude, that apart from degeneracies 
following from the single dot [that is limits d — ► 0(oo)] 
and possible accidental degeneracies respecting the selec- 
tion rule, double dot states are described by an effective 
Hamiltonian 

Hf = E?»-±Eac&. (29) 

Particularly, the energies of the two lowest states are 
given by this equation, with an exception for the state 
r^i in the region of small d where it is coupled to T] 1 
through l z and we have to describe it here by a 2 x 2 
effective Hamiltonian. 

An illustration of the l z influence on the spectrum is in 
Fig. [7| where the linear Dresselhaus spin-orbit contribu- 
tion to the energy for several states as a function of the 
interdot distance is shown. One can see at what interdot 
distances the l z operator causes the qualitative change 
between the double dot case (where the functions are 
characterized by a definite representation r, and the en- 
ergy contribution from the spin-orbit is a uniform shift) 
and the single dot case (where the functions are num- 
bered according to the orbital momentum and the spin- 
orbit contribution to the energy depends on a z l z ). This 
happens when E^ajj is comparable to the energy dif- 
ference of the nearly degenerate states. If the criterion 
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FIG. 7: Calculated corrections to selected lowest energy levels 
due to Hd- All states have spin a — +1. The graph demon- 
strates a transition between the symmetry group of the dou- 
ble dot Ho (states F) and that of single dots (states ^). The 
transition is induced by l z which by symmetry couples states 
Ti <-> F3 and Ti <-» IV Thus the anti-crossing mechanism will 
induce transition F^) <-> Ti ±F3 and F 2 (4) <-> F 2 ±T4. These 
linear combinations are equal to a single dot solution ^ n ,i,a in 
the case d — > and a combination a = &ni a ± a of 
functions with the same orbital momenta in the case d — > oo. 

for the coupling between the dots is the constant con- 
tribution, —a D Eo/2, to the energy, then the double dot 
region, as far as the spin-orbit coupling is concerned, is 
between 1 to 100 nm, that is up to 5 times of the confine- 
ment length of 20 nm. As an example, for the function 
r^ 1 the coupling in the effective Hamiltonian through l z 
to r^ 1 is equal to the unperturbed energy difference if 
a D ~ d 3 e~ d , giving d w 3, corresponding to the inter- 
dot distance of 6Zo- Due to the exponential, this result is 
insensitive to an- 

The Bychkov-Rashba term can be treated analogously. 
The effective Hamiltonian is Hff = [H - (a 2 BR /2)E (l + 
a z l z )]ij. The absence of a linear influence on the energy 
was based on the existence of H D P '. Since we found a 
case where -ff£>3 causes linear anti-crossing (see discus- 
sion to Fig. 0}, it follows that H^ can not exist for 
our double dot potential. However, if one approximates 
Ei — Ek ~ Ej — Ek in 1)25(1 . one can use if£f to sim- 
plify the mixed D-D3 correction. This, according to Fig. 
is the dominant spin-orbit correction for the tunnel- 
ing energy SE t . One gets an analogous expression as Eq. 
(|26|l . where the needed commutator is stated in Eq. JTSJ. 
Concluding, if we neglect the mixed D3-D3 term, we can 
write the spin-orbit contribution to the energy for the 
lowest orbital states to be (i — S, A) 

5Ef° = -E {a 2 D + a 2 BR )/2 + a D a D3 E l 2 Q (k 2 ) H . (30) 

One note to the eigenfunctions: The matrix elements of 
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the effective Hamiltonian are computed using the eigen- 
functions of Hq. But the functions that correspond to 
the solutions are transformed by the same unitary trans- 
formation that leads from Hq to H cfi . The sum rule can 
be used also here to express the influence of H\i n on the 
eigenf unctions of Hq. If HoTi = B,r,, the eigenf unctions 
corresponding to the effective Hamiltonian, Eq. I|26|l , are 



(31) 



Partial diagonalization of the effective Hamiltonian, to go 
from Eq. I|2fci|) to Eq. J2HJl, means we finish the unitary 
transformation completely and get Tf = Tf — H°^Tf for 
the eigenfunctions corresponding to the effective Hamil- 
tonian, Eq. 



the LCSDO approximation: 



E 



(o) _ 



E 



2E ( 1 + [1 - d V (l - 2 )/V^]e-^) 2 ( 1+92 ) 



S V 2 



1 + e -(*?) 2 (l+0 2 ) 

dt){l - 6 2 )[e-^' 2 - d7?Erfc(rfr?)] \ 
1 + e -(d v ) 2 (i+0 2 ) )■ 

(0) = 2Eo (I - [1 - dr,(l - e 2 )/0F]e-( d ") 2 ( 1 + 92 ) 

A ^2 



1 - e -(*?) 2 (l+0 2 ) 

dr)(l - 2 )[e-W /a/tt - tfy Erfcfefy)] 
1 _ e -(d»j) 2 (i+e 2 ) 



.(34) 



Here r\ — ?o/Zb = (1 — 2 ) 1 / 4 . In the limit <i — > oo, we can 
then deduce the tunneling energy in the leading order to 
be 



SB® = Eo-$=(l - 0^ e -d 2 (i+0 2 )/Vr- 



(35) 



C. Finite magnetic field. 

The presence of a magnetic field lowers the symmetry 
of the Hamiltonian without spin-orbit terms. The only 
nontrivial symmetry operator is the inversion / (see Tab. 
[nj. As a consequence the double dot states fall into two 
groups (representations of C 2 ): I\ and T 3 become 
(symmetric under I) and r 2 and T4 become Ta (anti- 
symmetric under /). Symmetrized functions <?™' i,CT now 
are 



n,l,a <Tr — d 1 7-1 ,Tfd 
9i =^n,l,* + A*„,J )<T , 



(32) 



where the irreducible states i — S and A, while the per- 
mutation coefficients D$ = —Da = 1. The shifted single- 
dot wave functions acquire a phase: 

V) = ± l od, y)e ±m ^Hl , (33) 

depending on which dot they are located. 

The double dot energy spectrum of Hq as a function 
of magnetic field is shown in Fig. [S] for the interdot dis- 
tance of 50 nm. Indicated are two crossings and one 
anti-crossing induced by magnetic field. The first cross- 
ing is between Ti and r 2 (notation from the B = case). 
These two states have opposite spins so they are not cou- 
pled and there is no level repulsion here. (We will see in 
the next section that spin-orbit coupling will induce anti- 
crossing in this case.) The second crossing is between r 2 
and T3, which behave differently under / and so they are 
not coupled by magnetic field. The actual anti-crossing 
is between T 2 and T4, which are both antisymmetric un- 
der /. This is an example of anti-crossing induced by 
magnetic field. 

In analogy with Eq. 11' I'D we derive analytical expres- 
sions for the energies of the lowest symmetric and anti- 
symmetric states in the presence of magnetic field using 



If 9 = 0, the above expressions reduce to Eq. (|22I23[1 . On 

B -b/2 e -B/B _ 



the other hand, if B 



o, then 5E t 

Magnetic field suppresses SE^ by suppressing overlap 
integrals (&- d \H\H! d ). There are three different effects 
that magnetic field introduces. First, the wave func- 
tions are squeezed by the confinement provided by the 
vector potential. The natural confinement length is 
lg = Zo ( 1 — 2 ) 1 / 4 , present in the exponential decay fac- 
tors. Second, the gauge phase produces factors (1 + 8 2 ) in 
the exponents of the scalar products (\& _d |\I fd ). Third, as 
B increases, the confinement potential Vc becomes less 
important compared to the confinement of the magnetic 
vector potential. This gives rise to the factor (1 — 2 ). 
Note that in the limit B — > 00, ; a is an eigenstate 
of Hq for any d. For reasons explained in last para- 
graph in the Sec. IIV Al Eqs. I)34|) are correct in the 

limit d — ► 00. As d —> 0, only E^ is correct. The limit 



value of E A °\d = 0)/E is 4/(1 + 6> 2 ), instead of the exact 
value of 2(2-61). 

At a finite magnetic field we have also a new term in 
the Hamiltonian, the Zeeman term. Since it commutes 
with Hq the only consequence of this term is a shift of 
the energy of the states by a value v/ibB according to 
their spin a. Therefore it introduces new crossings of the 
states with opposite spin. An example of this can be seen 
in Fig. where we plot energies of the four lowest states 
in the region where the Zeeman shift is comparable to 
the energy differences of the considered states. 



D. Effective spin-orbit Hamiltonian 

We now study the influence of spin-orbit coupling on 
the spectrum of double dots in a finite magnetic field. 
We will see that spin-orbit terms lead to new spin hot 
spots even at magnetic fields of the order of 1 T, and 
that linear spin-orbit terms will influence tunneling in 
the second order. 




FIG. 8: Computed energy spectrum of the double dot Hamil- 
tonian without Hso, as a function of magnetic field. The 
quantum dot separation is 50 nm (single dot confinement 
length is 20 nm). The energy levels are labeled according 
to the symmetry of the states at B — 0. Two crossings (one 
between T2 and T3, the other between Fi and F2) and one 
anti-crossing (between F2 and I^) are indicated. In the limit 
B — > oo the states merge to Landau levels. 

Although the presence of the Zeeman term complicates 
the analysis of the perturbation theory using operators 
H op , one can still apply the previously developed formal- 
ism if the Zeeman term is treated as a part of perturba- 
tion. (For a harmonic potential describing single dots, 
operators have been derived^ for the case of finite 
magnetic field, so that the Zeeman term can be included 
into Hq). Up to the second order in the perturbation 
couplings (being now a so and a z), there is no coupled 
Zeeman-spin-orbit term. This means that in the effec- 
tive Hamiltonians H cS that we already derived for the 
case of zero magnetic field, the Zeeman term appears as 
a shift of the energies on the diagonal without bringing 
any new couplings (non-diagonal terms). But an impor- 
tant consequence is that the shift can change the number 
of states we have to include into the basis where the ef- 
fective Hamiltonian acts, because their energy difference 
to the considered state is comparable to the spin-orbit 
coupling. 

First, in analogy with Eq. (|29(l . if the energy of a state 
is far enough from others, we can consider the basis to 
consist of one term only and the spin-orbit correction to 
the energy of state \i) is 

+[H°*,H m ], (36) 

where the averaging means the expectation value in the 
state \i) and a is the spin of the state. Since the pres- 
ence of magnetic field lowers the symmetry, the last com- 
mutator, [Eq. I|15|) ]. can not be simplified according to 



FIG. 9: Calculated energies of the four lowest states of Hamil- 
tonian Ho + Hz at B = 1 T. The vertical dashed lines indi- 
cated doublets in the effective Hamiltonian. The inset dis- 
plays the anti-crossing at 89 nm due to Hd- Dashed lines are 
energies of Ho + Hz, solid lines of Ho + Hz + Hd- 

the symmetry as was the case before in Eq. (|3U|1 . and, 
more important, we no longer have L z — 0. As a result, 
there are now corrections to the tunneling that are of the 
second order in the linear spin-orbit couplings. These 
corrections depend on a_ = a 2 D — a 2 BR . 

Second, we look how the energies of the four lowest 
states are changed, using again the example of the linear 
Dresselhaus term. They are plotted in Fig. 03 Here in the 
main figure one can see the shift caused by the Zeeman 
term and the anti-crossing induced by the spin-orbit cou- 
pling is magnified in the inset. The anti-crossing states 
are ri and v\. In the case of zero magnetic field we 
described each of the four basis states by Eq. I|29[) . Now, 
in principle, we have to describe them by a 4 x 4 effective 
Hamiltonian Eq. (|26|l . Due to symmetry we can simplify 
this Hamiltonian into two 2x2 Hamiltonians, i?i ff , ff| ff , 
acting in the bases r^,r^ and rj^r^ respectively. The 
four components of the effective Hamiltonian matrix are 

2 

(H cS ) n = E^ } -E ^-(l-a(L z ) u )-a^ B B-R n 

2 

(H cS ) 22 = Ef -E ^-{l+a(L z ) 22 )+afi B B + R n 

(H eS )i2 = (H cS )l = {H D ) U (37) 

where a = 1 for iJJ and a — — 1 for if| , while indices 
1,2 denote the first and the second term in the corre- 
sponding basis. Comparing to the case of zero magnetic 
field the Zeeman term increases the difference of the di- 
agonal elements in i?° ff and decreases them in ■ The 
ground and the fourth excited states which are described 
by Hf s stay isolated, and we can do the perturbative 
diagonalization to get rid of the off-diagonals. The en- 
ergy of the two states is then up to the second order in 
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FIG. 10: Calculated corrections to the energies of the four 
lowest states due to the linear Dresselhaus term Hd , at B = 1 
T. Solid lines are numerical data, dashed lines are analytical 
expressions computed by Eq. 136H using the LCSDO approx- 
imation for the states. 



the spin-orbit coupling accurately described by Eq. (J3SJ). 
Concerning the states Tg and T\, there is a region in 
the intcrdot distance of a few nanometers, where the two 
states must be described by the two dimensional iff to 
account for the anti-crossing, which is caused by the the 
matrix element (T^HdIT^). LCSDO gives for this cl- 
ement a result correct only in the order of magnitude. 
This is because even in the limit d — > oo this matrix el- 
ement is of the same order as the neglected coefficients 
c(n,l) in the LCSDO approximation, Eq. (|20|l . 

The spin-orbit corrections to the energies from Ho for 
the four lowest states as functions of the interdot dis- 
tance are in Fig. ^3 Also shown are analytical values 
computed by Eq. (|36|) . that is, ignoring anti-crossing. 
The scale implies that the corrections are of the second 
order in ao and for the states Tg and T^ A are enhanced in 
the anti-crossing region. The incorrect values at the limit 
d — > for functions T a A are because of reasons explained 
in the last paragraphs in Sec. IIV Al 



E. Spin-orbit corrections to the effective g-factor 
and tunneling frequency 



We next analyze spin-orbit corrections to the g-factor, 
5g = [SE(Tg) — 8E(y\)]/ [IbB, that characterizes the 
energy cost of a spin flip in the ground state, or the fre- 
quency of a spin precession. Another kind of oscillation is 
electron tunneling, when electron oscillates between the 
left and the right dot without changing its spin. The fre- 
quency of this oscillation, SE t /2h, is given by the energy 
difference SE t = E(T' A ) — E(T^ S ), Corrections to this en- 
ergy difference induced by the spin-orbit interaction are 
denoted in this paper as 8Ef° . 

First, we take a look at the spin-orbit corrections to 
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FIG. 11: Calculated spin-orbit corrections (relative to the 
conduction band value g*) to the effective g- factor, Sg = 
[SE(Tg) — SE(Tg)]/ hbB, as functions of magnetic field. The 
distance between the dots is 50 nm. Solid lines are numerical 
data, dashed lines are analytical values computed using Eq. 
13011 . Contributions come from linear spin-orbit terms (a, b), 
the mixed Dresselhaus correction from Ho and Hd3 (c) and 
the cubic Dresselhaus Hds correction (d). 



the ^-factor. Contributions in the second order of the 
spin-orbit couplings are shown in Fig. ^2 as functions 
of magnetic field at a constant interdot distance. The 
spin-orbit contribution to the g-factor in the double dot 
case has the same qualitative dependence on the mag- 
netic field as in the single dot case (see discussion to Fig. 
121 • However, at finite interdot distances, there is an en- 
hancing effect on the spin-orbit contributions. This can 
be seen in Fig. ^3 where at a certain magnetic field, the 
spin-orbit contribution to the g-factor is enhanced for a 
finite d compared to the case of d = 0(oo). We found 
numerically, that the enhancement can be up to 50% of 
the value of the correction in d = at magnetic fields of 
the order of 1 T. 

At the anti-crossing the spin-orbit contributions show 
cusps. At magnetic fields bellow the anti-crossing, the 
dominant spin-orbit contribution is D-D which reduces 
the conduction band g-factor by several percent. Con- 
tributions D-D3 and BR-BR are one order of magni- 
tude smaller. Using Eq. (|36() . that is ignoring the anti- 
crossing, we get for the contribution from the linear spin- 
orbit terms 



Sg(lin - lin) = — -^a^L z , 
fi B B 

where, in the limit d — > oo, 

Ll^eil + idhfe-^ 1 ^]. 



(38) 



(39) 



From Fig. ^2 one can see that the analytical result agrees 
with numerics. 

Finally, we look at the tunneling energy in the pres- 
ence of both magnetic field and spin-orbit couplings. The 
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FIG. 12: Calculated spin-orbit corrections to the tunneling 
energy SEt as a function of magnetic field. The interdot dis- 
tance is 50 nm. Solid lines are numerical data, dashed lines 
are analytical expressions computed by Eq. 13611 . 

spin-orbit corrections, as a function of magnetic field, 
are shown in Fig. ^] At zero magnetic field there is 
no contribution from the linear terms (result of the sec- 
tion I1VB|) and the dominant contribution is D-D3. Sim- 
ilarly to 8E[ , the corrections decay exponentially with 
increasing magnetic field. Anti-crossing strongly influ- 
ences the tunneling energy. Using LCSDO for d — ■> oo we 
obtain in the leading order 

6E?» = -£ a (2) #(^) 2 e- (d " )2(1+92) - (40) 

This analytical formula underestimates the corrections 
from the linear spin-orbit terms by a factor of ~ 3. Nev- 
ertheless the analytical expression for D-D3 is reason- 
ably good. In the magnetic field below anti-crossing, the 
relative change of the tunneling energy stemming from 
the spin-orbit terms is of order 10~ 3 . 

F. Tunneling Hamiltonian 

We now use our results to describe the influence of the 
spin-orbit interaction on the lower part of the spectrum. 
We restrict our Hilbert space on the four lowest states 
•^jS(A)' t ne eigenstates of the total double dot Hamilto- 
nian. Because of the transformation Eq. I|31|l these four 
states have neither definite symmetry with respect to in- 
version I, nor a definite spin in z-direction. fn this section 
we will denote them as spin 'up' and spin 'down' states. 
For description of a transport through the double dot it 
is useful to define the following left and right localized 
functions 

L°{R°) = ±={T S ±T<> A ), (41) 

where plus and minus holds for L and R, respectively. In 
the limit d — > oo these functions converge to single dot 



solutions localized in the left and right dot. 

The effective Hamiltonian of our system in the second 
quantization formalism is 

H= E a (n LrT + n Ra )~{SE^/2)(a[ a a Ra + a Ra a La ), 

(42) 

where E a = {E% + E%)/2, 6E? = E A — Eg, while a*, 
and a are creation and annihilation operators, and n — 
a) a. We can get both localized and spin-pure states if 
we diagonalize o~ z in a chosen basis. For example, taking 
and we get L purcT ~ (£t + oL l ) and L pure ^ ~ 
{L* — o + L^), up to normalization (f — |o| 2 /2). That is, 
the left pure spin state is a linear superposition of both 
left states with spin 'up' and 'down'. The coefficient o is 
proportional to esso- 
in the following we are interested in the time evolution 
of localized states given by Hamiltonian Eq. 1|42|) . First 
we note, that due to the non-diagonal terms, the electron 
which is in a localized state will tunnel into the other 
localized state after the tunneling time t% nn = h/5E^, 
resulting in coherent oscillations. For our parameters 
ttun ~ f P s - In the Hamiltonian there is no mixing be- 
tween spin 'up' and 'down' states. However, there will be 
mixing (or spin-flip) if we work with localized pure spin 
states. Electron being originally in £P ure T will, after the 
tunneling time t\ un , contain i?P uro ^ with the probability 
amplitude 

c = wk(5E\ - 5eJ ) /25Ej , (43) 

assuming that the difference in SE t for different spins is 
much smaller that 8E t itself. 

In the case of zero magnetic field, because of Kramer's 
degeneracy, the tunneling frequencies are the same for 
both spin orientations. Then whatever is the initial com- 
bination of spin 'up' and 'down' (let it be, for example, 
£puret ) ^ cl ur ing the time evolution (oscillations) there will 
be no relative change in the coefficients in this linear com- 
bination. Therefore spin-orbit coupling does not lead to 
spin- flipping and c = in Eq. (|43|l . 

In a finite magnetic field, the tunneling frequency for 
spin 'up' and 'down' are different. The difference is 

(2) 

caused only by spin-orbit terms, and is of order a_ = 
cvq — cCqj^. Equation l|3()|) shows, that the contribution to 
5E\ from the linear spin-orbit terms is opposite that of 
5E\ and therefore their difference is twice the expression 
in Eq. (|40|l . We conclude, that spin-flip during tun- 
neling induced by spin-orbit coupling is proportional to 
the third power in spin-orbit couplings and depends lin- 
early on the magnetic field if the magnetic field is small 

The different tunneling frequency can be exploited for 
separation of different spin states in a homogeneous mag- 
netic field. Starting with some combination of 'up' and 
'down' states localized in one dot, after time t sop — 
h/(SEj-SE^) = U u JE}/{5E}-8E l t ) the part with spin 
'up' will be localized in the left, and the spin 'down' will 
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be in the right dot. From Fig. ^]one can see that one 
needs about 10 3 coherent oscillations to get the spatial 
spin separation. Therefore the decoherence time must 
be longer in order to observe this effect experimentally. 
We note that the separated states will not be pure spin 
states, but will contain a small (linearly proportional to 
aso) admixture of opposite pure spin states. 

V. CONCLUSIONS 

We have performed numerically exact calculation of 
the spectrum of a single electron localized by a confin- 
ing potential in single and double GaAs quantum dots. 
We have studied the influence of the spin-orbit terms, 
namely the Bychkov-Rashba and the linear and cubic 
Dresselhaus terms, on the energy spectrum. In the sin- 
gle dot case we have elaborated on previous results and 
shown that the spin-orbit interaction has three principal 
effects on the spectrum: first, the interaction shifts the 
energy by a value proportional to the second order in 
the spin-orbit couplings, second, it lifts the degeneracy 
at zero magnetic field, and, third, the Bychkov-Rashba 
term gives rise to spin hot spots at finite magnetic fields. 

In the double dot case we have addressed the symme- 
tries of the Hamiltonian. For zero magnetic field with- 
out spin-orbit terms we have constructed the correlation 
diagram, between single and double dot states, of the 
spectrum. We have used properly symmetrized linear 
combination of shifted single dot solutions as an approx- 
imation for a double dot solution and found that for the 
four lowest states it gives a good approximation for the 
energy. As for the contributions to the energy from the 
linear spin-orbit terms, we have found that in zero mag- 
netic field a typical feature of a double dot is a uniform 
shift of the energy proportional to the second order in the 
coupling strengths without any dependence on the inter- 
dot distance. This is true also if the potential has lower 
(or none) symmetry (for example biased dots). There- 
fore, in zero magnetic field, there is no influence on the 
tunneling frequency up to the second order in the lin- 
ear spin-orbit couplings and the dominant contribution 
comes from the mixed linear and cubic Dresselhaus sec- 
ond order term. We found, that spin hot spots in zero 
magnetic field exist in the double dot, but are solely due 
to the cubic Dresselhaus term. This means also, that 
for our potential, for the cubic Dresselhaus term there 
can not exist a unitary transformation to eliminate its 
contribution in the first order. 



The effective (/-factor, on the other hand, is influenced 
by the second order linear spin-orbit couplings even at 
B rsj 0, so the dominant contribution here is the linear 
Dresselhaus term for GaAs. In finite magnetic fields the 
uniform shift does not hold any more and there is a con- 
tribution to the tunneling frequency in the second order 
of the linear spin-orbit couplings. We have derived an 
effective Hamiltonian, using Lowdin's perturbation the- 
ory, with which analytical results up to the second order 
in perturbations (Zeeman and spin-orbit terms with the 
exception of cubic Dresselhaus-cubic Dresselhaus contri- 
bution) can be obtained provided one has exact solutions 
of the double dot Hamiltonian without Zeeman and spin- 
orbit terms. From this effective Hamiltonian we have 
derived the uniform shift in zero magnetic field. In a 
finite magnetic field we used linear combinations of sin- 
gle dot solutions to obtain analytical expressions for the 
spin-orbit contributions to the energy for the four low- 
est states. We have analyzed them as functions of the 
interdot distance and magnetic field and compared them 
with exact numerical values. The spin-orbit relative con- 
tribution to the g-factor and the tunneling frequency is 
of the order of ~ 1CP 2 and ~ 1CU 3 , respectively. Due to 
the degeneracy of the energy spectrum at large interdot 
distance the spin hot spots exist also at smaller magnetic 
fields compared to the single dot case. 

As an application of our results we have constructed 
an effective Hamiltonian acting in a restricted Hilbcrt 
space of four states-electron with spin up and down 
(these arc effective spins in the presence of spin-orbit cou- 
pling) localized on either dot. Effectively, there is only 
spin-conserving tunneling between the localized states, 
no spin-flip tunneling. In zero magnetic field the spin- 
orbit interaction does not significantly influence the tun- 
neling frequency, nor it implies spin-flip tunneling. In 
finite magnetic fields the tunneling frequency is spin de- 
pendent, the difference being of second order in linear 
spin-orbit terms. This leads to a spin flip amplitude pro- 
portional to the third power in spin-orbit couplings (it is 
linear in magnetic field). We propose to use this differ- 
ence of the tunnelings to spatially separate electron spin 
in homogeneous magnetic field. 
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